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Abstract 

A new method of calculating pairing correlations in coordinate space with 
finite range interactions is presented. In the Hartree-Fock-Bogoliubov (HFB) 
approach the mean field part is derived from a Skyrme-type force whereas 
the pairing field is constructed with a Gogny force. An iterative scheme is 
used for solving the integro-differential HFB equations via the introduction 
of a local equivalent potential. The method is illustrated on the case of the 
nucleus 18 C. It is shown that the results are insensitive to the cut off value in 
the quasiparticle spectrum if this value is above 100 MeV. 


The treatment of pairing correlations is very important for the description of the prop¬ 
erties of weakly bound nuclei situated close to the drip lines. The Hartree-Fock-Bogoliubov 
(HFB) method [|T| is the the commonly adopted approach for treating self-consistently both 
the mean held contributions and the pairing correlations. Various effective interactions are 
available and can be used in the mean held and pairing channels of the HFB equations. 
They are of two types, namely zero-range forces |2j-[|| and finite range forces whose typical 
representatives are the Gogny forces [0. 

Zero-range forces are widely used because the self-consistent equations can be conve¬ 
niently solved in coordinate space. In loosely bound systems such as unstable nuclei, contin¬ 
uum effects in the pairing channel become important since their Fermi energies are typically 
close to zero and the pairing correlations can easily populate the continuum. Then, working 
in coordinate space is an advantage because the continuum effects can be accurately treated 
On the other hand, a zero-range pairing interaction has the well-known pathology of 
producing diverging contributions if no cut oh is imposed on the quasiparticle space. This 
cut oh must be an inherent part of the phenomenological zero-range interaction [10], but it 
is not clear which cut oh value must be adopted for a given interaction and a given nucleus. 
We note that some regularization scheme has been proposed recently |TT[ for dealing with 
the question of HFB equations with zero-range pairing interactions. 

Finite range pairing interactions of course do not require in principle a truncation of the 
quasiparticle space, even though in practical calculations the summations are restricted to 
quasiparticle energies below some cut oh energy E c o . If E co is large enough the results 
no longer depend on its precise value, in contrast to the case of zero-range interactions. 
It is desirable therefore to have a method which combines the advantages of solving HFB 
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equations in coordinate space and of describing pairing correlations with a finite range 
interaction. It is the purpose of this paper to propose such a method and to illustrate it by 
comparing results obtained with Skyrme and Gogny interactions. 

We treat the continuum states exactly as explained in Ref. || where Skyrme-type inter¬ 
actions were used both in the mean field and pairing channels. Here, we keep a Skyrme force 
for the Hartree-Fock (HF) mean held but the pairing interaction is taken as a Gogny force. 
The calculations are done for the neutron-rich nucleus 18 C with the Skyrme interaction SLy4 
IP in the mean held channel and the Gogny interaction D1S [[E| in the pairing channel. 

We would like to stress here that with these calculations we do not want to make any 
predictions on the properties of the chosen nucleus: to do this a ht of the parameters of the 
Skyrme force to be used with the Gogny pairing force must be done. However, the purpose 
of this work is not to make predictions on the properties of the nucleus 18 C but to present 
and show the advantages of a method of solving HFB equations where a zero-range and a 
finite range interactions are combined together. 

For clarity we recall the general form of the coupled integro-differential HFB equations 
in coordinate representation p[: 


fdW ( h{ra,r'a' \ ( 0i(£,rV)\ 

a ' \ h(ra,r'a') —h(ra,r'a'J \4> 2 {E,r'a') J 

( 1 ) 

f E + E f 0 \ f (j)i(E,ra) \ 

V 0 E-E f ) \(^(E,T(7)) ’ 


where E F is the Fermi energy, (j>i and (j ) 2 are the upper and lower components of the HFB 
quasiparticle wave functions, respectively. In Eq.(|l|) the HF operator is a sum of kinetic and 
mean field components: 


hirer, rV) = T(r, r')<W + r(rcr, rV) , (2) 

r(rcr,rV) = f d 3 r 1 d 3 r 2 V(ra,r 1 a 1 -r'a , ,r 2 a 2 )p(r 2 a 2 ,r 1 ai) . (3) 

When the HF interaction V is a Skyrme force, the mean field T is just a functional of the 
local particle density p( r) and of its derivatives. Then, the HF operator h is a differential 
operator, which is a familiar property of Skyrme-HF models. On the other hand, the pairing 
interaction V pa i r is a local force of finite range, and the pairing operator h(ra, r 'a') remains 
a fully non-local kernel obtained by folding V pair with the pairing density p\ 

h(ra, rV) = / d 3 r 1 d 3 r 2 ^ 2a'a' 2 V pair {ra, r' - </; ricq, r 2 - a 2 )p{r 1 a l) r 2 cr 2 ) . (4) 

This is the new feature as compared to the equations discussed and solved in Refs. [!!• 
We choose to work with positive energy quasiparticle states, and the particle and pairing 
densities are expressed as: 

p(rcr, rV) = ^(F’n, rcr)(j)* 2 {E n , r'cr) + f dE</) 2 (E,ra)(j)l(E,r'a) , (5) 

0 <E n <-E F j ~ Ef 
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rEc.o. 


p(rcr, r'a) = — Y 4>2(E n , ra)(j>l(E n , r'a) - / dEfciE, ra)4>\(E, r'cr) 


0<^ n <--^F 


—Ep 


( 6 ) 


In Eqs. (|5|) and © the summations are over the discrete states that are situated in the 
region of the spectrum 0 < E n < —Ep whereas the integrals run over the continuum part 
of the spectrum up to a chosen cut off. 

We now take for V pa i r a Gogny force which contains a sum of two gaussians, a zero-range 
density-dependent part and a zero-range spin-orbit part. Within the parametrisation D1S 
m that we adopt in this work the zero-range density-dependent part does not contribute. 
Let us explain the contribution of the gaussian terms. That of the spin-orbit part is also 
included in the calculations presented below although its effect is quite small. The finite 
range part of the interaction is: 


Vpair( r i - r 2 ) = Y ( w ° + B aP<r - H a P T - M a P (T P T )e 

a= 1,2 


— I r l — r 2 I ~ 

m2 


( 7 ) 


where XV Q , B a , H a , M a and p a are parameters, P CT and P r are the spin and isospin exchange 


operators, respectively. Then, the pairing operator (ffl) becomes [10 


h(ra, r'a ) = Y e ^ [(W a - H a )p(ra , r'cr) - (B a - M a )p{r'a , rcr)] . (8) 

a=l,2 

In the following we will need the multipole expansions of the gaussian form factors: 

=EF£(r,r')Y LM (f)Yl M (r’) . (9) 

LM 

We restrict ourselves to the case of spherical symmetry, for simplicity. Then, the general 
set of equations © can be solved for each partial wave ( l,j ) separately. We define the radial 
parts of the quasiparticle wave functions &i(E,ra) (i=l,2) by: 


<*i (Eljm,ra) = ‘EVt — Y lm (f)(lm,^cy\jm)x(a) , 


( 10 ) 


where \(cr) is a spinor corresponding to a spin projection a. With the help of the definition 
(|TC|) and the expansion ([|) it is straightforward to obtain the multipole decomposition of 
the pairing field h{ r, r') = h(ra, r'a): 


where 


Hr,r')= Y h Ll (r,r')Y LlMl (r)Y£ lMl (f') 

LxMx 


h Ll (r, r') = Ea= 1)2 El Fftr, r')^±i £„y(2j + 1) 

(]{ — W \ u 2 (nlj,r) m(nZj,r') _ g \ u 2 (nlj,r') ui(nlj,r) 

\ Gi J rr*f \ GL Q! / ry* 

( L l LA 2 

l 0 0 0 j ' 


( 11 ) 


( 12 ) 
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The summations over n in Eq. (|T2D become integrals over the energy 
states. 

For each partial wave ( Ij) one has to solve a system of two coupled 
equations whose general structure is: 


for the continuum 
integro-differential 


hui(r) + J h(r,r')u 2 (r')r' 2 dr' — (E + E F )u\{r ) , 

J h(r, r')ui(r')r ,2 dr' — hu 2 {r) — (E — E F )u 2 [r) . (13) 


In the context of HF equations with finite range interactions it has been shown by Vautherin 
and Veneroni |T^j that one can transform the HF integro-differential equation into a purely 
differential equation by introducing a so-called trivially equivalent local potential. Here, we 
generalize this method to the system of equations ([13]) by defining local equivalent potentials 
Ui(r) in the following way: 


h(r, r')ui{r')r ' 2 dr ' = 


Ui(r) 


h(r, r')ui(r')r' 2 dr')ui(r) 


= Ui(r)ui(r) , i = 1,2 . 


(14) 


Then, Eqs.(|T3|) become formally a system of two coupled differential equations where the 
potentials depend on the solutions and therefore they must be solved iteratively. This is not 
a major problem since the self-consistency requirement already leads to an iterative scheme. 

An additional difficulty comes from the fact that the local potentials U t (r) have poles 
at the nodes of the wave functions Ui(r). In Ref. [II] a very simple and efficient method 


was proposed to overcome this problem, based on the linearization of the local equivalent 
potential around the poles. Another equivalent potential LA(e, r) is introduced; it is equal to 
Ui{r) everywhere except inside the intervals [r 0 — e, r 0 + e] where r 0 denotes a pole of Ui(r). 
Inside these intervals Lj(e, r) is chosen as a segment which joins the values of Lj( r o — e ) and 
Ui(r 0 + e). The approximation is good if e is small enough not to wash out the shape of the 
potential. Thus, Ui(r ) is replaced by the following potential: 


U^ n+1) (e n+1 ,r) = R £n+1 ( [ r ,2 dr'h{r,r')ui(r') {n) dr'\ , (15) 

l u\ '(r) J J 

where R en+1 indicates the linear interpolation procedure described above and n counts the 
iterations by which the HFB equations are solved. The parameter e depends on the iteration 
and it is chosen so that lim e n = 0. At each iteration of the HFB scheme we evaluate the 
equivalent potentials of all quasiparticle states by using the wave functions of the previous 
iteration and we repeat this procedure until convergence. 

We compare now the HFB results for 18 C obtained by using in the pairing channel 
either the finite range Gogny interaction D1S or a zero-range interaction. In both cases the 
Skyrme interaction SLy4 is used to construct the mean field. The quasiparticle continuum 
is fully treated (no box boundary conditions) as described in Ref. [|9|. As it was said before, 
the problem with the use of a zero-range pairing interaction is the divergence of the pairing 
correlations when one increases the cut off energy F c . 0 .. This is illustrated in Fig.l where the 
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particle (top) and pairing (bottom) densities of neutrons in 18 C are plotted. The zero-range 
pairing interaction in this case is 0 : 



r') 


V 0 5(r - r') 



(16) 


where p(r) is the particle density, and the values of Vo and 7 are fixed from the parameters 
of SLy4 (Vo=-2488.9 MeV fm 3 , 7 = 1 / 6 , with the choice pc = 0.133 fin -3 ). This particular 
choice of the parameters does not influence significantly the convergence properties at high 
energies cut-offs discussed here. 

In Fig.l different densities corresponding to different values of E c _ 0 _ are plotted. We can 
observe that, while the particle density is almost stable with respect to the cut off energy 
the pairing density increases with E CX) . as it is expected. This indicates that the mean field 
properties are not much affected by the enlarging of the continuum phase space while the 
pairing properties are very sensitive to it. If one moves the cut off towards higher values the 
pairing density continues to increase and it never converges to a stable result. This problem 
is eliminated when a finite range interaction is used to construct the pairing field. 

We show in the upper part of Fig.2 the total energies of 18 C calculated with the two 
pairing interactions and for cut off values ranging from 70 MeV up to 120 MeV. We observe 
that the system is more and more bound when the cutoff increases; this is due to the fact that 
pairing correlations become more and more important. While in the case of the zero-range 
interaction the total energy continues to decrease in the chosen interval of cut off values, in 
the case of the Gogny interaction the total energy converges and reaches a stable value equal 
to -130.65 MeV at a cut off of 100 MeV. We can equivalently observe this convergence of 
the energy by studying the trend of the pairing correlation energies which are shown in the 
lower part of Fig. 2 for the same cut off values. The pairing correlation energy is defined as 
follows: 


E P = E(HF) - E(HFB) , (17) 

where E indicates the total energy of the nucleus. The quantity Ep is the difference in the 
total binding energies between the HF and HFB calculations; thus, it gives an estimation of 
the amount of pairing correlations. In the figure one can observe that these energies always 
increase for the zero-range interaction indicating the increasing of the amount of pairing 
correlations while they reach a stable value equal to 16.04 MeV for the Gogny interaction. 
It is easy to understand why the results with the latter interaction become stable for cut off 
energies around 100 MeV or beyond. The shortest range of the two gaussian form factors is 
0.7 fm, which corresponds to 2.86 fm -1 in momentum space, i.e., a kinetic energy of about 
160 MeV. The depth of the mean field potential being about 40 MeV one can estimate that 
this kinetic energy corresponds to a quasiparticle energy around 120 MeV. 

Another quantity that one can study as a function of the cut off energy is the mean 
square radius of the pairing density. Indeed, the mean square radius of the particle density 
must be stable against E c , 0 _ because the u\{r) functions entering Eq.([|) are negligible for 
large quasiparticle energies. In contrast, the ui(r)u 2 (r) factors of Eq.(| 6 |) do not decrease so 
fast with increasing quasiparticle energies. Let us define: 

( r a ) = j r 2 drp{r)r a . (18) 
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In Fig. 3 we show (r 2 ) for the two pairing interactions. Again, while (?’ 2 ) always increases in 
the case of the zero-range interaction it reaches a stable value of about 34.5 fm 2 at a cut off 
of 95 MeV in the case of the Gogny interaction. This convergence indicates that the pairing 
density does not change when the cut off is moved beyond the value of 95 MeV. 

In this work we have presented a new method for solving HFB equations in coordinate 
space with finite range pairing interactions. This may be useful in systems where the chem¬ 
ical potential is close to zero and an accurate treatment of the quasiparticle continuum is 
required. As an illustration of the method we have solved the HFB equations for a neutron- 
rich nucleus, using the finite range force D1S as pairing interaction. Because our main 
purpose here is to discuss the respective behaviour of zero-range and finite range interac¬ 
tions in the pairing channel, we have kept in the present application the Skyrme-type force 
SLy4 to generate the HF mean field. A further step toward full self-consistency with the 
same finite range force for calculating the mean field and the pairing field can be made, using 
the same technique of local equivalent potential also for the mean field as it was already 

0 . 


done in the HF context 


Work in this direction is in progress. 
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FIGURES 


FIG. 1. Neutron particle (top) and pairing (bottom) densities of 18 C obtained with the 
zero-range pairing interaction, Eq. (|l6|), for different values of the cut off. 


FIG. 2. Upper part: Total energies of 18 C obtained with zero-range and Gogny pairing inter¬ 
actions for different cut off values. Lower part: Pairing correlation energies, Eq. (0) , obtained 
with zero-range (circles) and Gogny (stars) pairing interactions for different cut off values for 18 C. 


FIG. 3. Mean square radii of pairing densities, Eq. (18), obtained with zero-range (circles) and 
Gogny (stars) pairing interactions for different cut off values for 18 C. 
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